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We present the result of a search for charge asymmetry in B ± — > J/tpK ± decays using 772 x 10 6 
BB pairs collected at the T(45) resonance by the Belle detector at the KEKB asymmetric-energy 
e + e~ collider. The CP-violating charge asymmetry is measured to be Acp(B + — > J/tpK + ) = 
[-0.76 ± 0.50 (stat) ± 0.22 (syst)]%. 

PACS numbers: 11.30.Er, 13.25.Hw 



Violation of CP symmetry in the Standard Model 
(SM) has been well established. Interest has now shifted 
to the search for new sources of CP violation due to 
physics beyond the SM, since the CP violation content of 
the SM does not explain the matter-antimatter asymme- 
try of the Universe [l[. In the SM, CP- violating phenom- 
ena in the quark sector are described by the Kobayashi- 
Maskawa theory Q , in which a single irreducible complex 
phase gives rise to all CP- violating asymmetries. The 
decay B + — » J/ipK + is mediated by a color-suppressed 
b — > ccs transition, where the dominant b — > c tree- 
level amplitude and the b — > s penguin amplitude have 
a small relative complex phase, arg [— V cs V* b /V ts V t l] @ 
in the SM Q. The charge asymmetry in the decay 
B + -> J/ipK + is defined as 

A C p{B+ ->■ J/^K+) 

_ B{B- -» J/i/jK-) - B(B+ -> J/i/jK+) 
~ B{B- -> J/4>K-) + B{B+ -> J/i/jK+) ' 

where B denotes the branching fraction. Direct CP vi- 
olation would appear as a nonzero Acp(B + J/tbK + ) 
value and is predicted to be quite small, ~ 0.3% [4j, in 
the SM. Some new physics models predict enhanced val- 
ues of this asymmetry For example, a model with 
an extra U(l)' gauge boson [H[ and another model with 
an extra coupling to the charged Higgs boson [|| predict 
asymmetries of 0(1%) and O(10%), respectively. 

The B + — > J/ipK + decay mode has low backgrounds 
and a large branching fraction, when compared to that 
of other charmonium decay modes. These properties en- 
able a precise Acp(B + — > J/ipK + ) measurement, which 



provides an important test of various new physics models 
and constrains their parameter spaces. 

The current world average for Acp(B + — > J/ipK + ) is 
(0.9 ± 0.8)% which is dominated by the DO result, 
(0.75 ± 0.61 ± 0.30)% 0, while the most precise result 
from an e + e~ collider experiment is the BABAR result 
of (3.0 ± 1.4 ± 1.9)% 0. 

In this paper, we report a measurement of Acp(B + — > 
J /ipK + ) using a 711 fb _1 data set that contains 772 x 
10 6 BB pairs. The B + — s- J/ifjK + decay is reconstructed 
in the J/tp — > (I = e or fx) channels. For a precise 

measurement of the charge asymmetry in this decay, the 
asymmetry of charged kaon detection efficiencies must 
be carefully studied and corrected for. The asymmetry 
in detection efficiency arises due to the asymmetric geom- 
etry of the detector, different interaction rates of kaons 
in the detector material, and differences in kaon identifi- 
cation efficiencies for K + and K~ . 

KEKB is an asymmetric electron-positron storage ring 
that collides 8.0 GeV electrons with 3.5 GeV positrons 
at the T(4S) resonance (center-of-mass [cm.] energy 
Vs = 2£ b cam = 10.58 GeV). The T(4S) resonance is 
boosted by #7 = 0.425 0. The Belle detector is a 
large-solid-angle magnetic spectrometer that consists of 
a silicon vertex detector (SVD), a 50-layer central drift 
chamber (CDC), an array of aerogel threshold Chcrenkov 
counters (ACC), a barrel-like arrangement of time-of- 
ffight scintillation counters (TOF), and an electromag- 
netic calorimeter (ECL) comprised Csl (Tl) crystals lo- 
cated inside a superconducting solenoid coil that provides 
a 1.5 T magnetic field. An iron flux- return located out- 
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side of the coil is instrumented to detect K\ mesons and 
to identify muons (KLM). The detector is described in 
detail elsewhere [ill EH ■ 

Hadronic events are initially selected by requiring at 
least three reconstructed charged tracks, a total recon- 
structed ECL energy in the cm. in the range O.ly/s — 
0.8-y/s, at least one large-angle cluster in the ECL, a total 
visible energy - calculated from all charged tracks and 
isolated neutral showers - greater than 0.2y/s, an ab- 
solute value of the z component of the cm. momentum 
less than 0.5-^/i, and a reconstructed primary vertex that 
is consistent with the known location of the interaction 
point. To suppress two-jet non-T(4S') background rela- 
tive to BB events, we require R2 < 0.5, where R2 is the 
ratio of the second to zeroth Fox- Wolfram moments [l3| . 
To remove charged particle tracks that are poorly mea- 
sured or do not come from the interaction region, we 
require dz < 5 cm for all tracks, where dz is the absolute 
value of the coordinate along the beam direction at the 
point on the track closest to the origin. 

The J I tp meson is reconstructed from one tightly and 
one loosely identified lepton. For muon tracks, the 
tight identification criterion is based on track penetra- 
tion depth and hit scatter in the KLM system, while 
the loose identification criterion requires only that the 
tracks have an energy deposit in the ECL that is consis- 
tent with that of a minimum ionizing particle. Electron 
tracks are tightly identified by a requirement on a com- 
bination of dE/dx from the CDC, E/p (E is the energy 
deposit in the ECL and the p is momentum measured 
by the SVD and the CDC), and shower shape in the 
ECL. For loose identification, either dE/dx or E/p is 
required to be consistent with the electron hypothesis. 
We correct for final-state radiation or bremsstrahlung in 
the inner parts of the detector by including in the e + e~ 
invariant mass calculation the four-momentum of every 
photon detected within 0.05 rad of the original electron 
or positron direction. Since small residual radiative tails 
still remain, we use asymmetric invariant mass require- 
ments, -60 MeV/c 2 < M^- -m J/4> < 36 MeV/c 2 and 
-150 MeV/c 2 < A'/ e+e - (7 ) - m J/4 , < 36 MeV/c 2 , for the 
/i + p~ and e + e~ pairs, respectively. 

The combined information from the CDC , TOF, and 
ACC is used to form a K — tt likelihood ratio, TZk = 
£k/(£k+£tt), where Ck(tt) is the likelihood of the kaon 
(pion) hypothesis. We require TZk > 0.6 for kaon can- 
didates, which is approximately 80% efficient for kaons, 
while giving a misidentification probability of below 10% 
for pions. 

Charged B mesons are reconstructed by combining a 
J/ip candidate with a charged kaon candidate. The en- 
ergy difference, AE = -Ecand — ^boam, and the beam- 
energy constrained mass, Mb c 
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used to separate signal from background {E ca _ n d and p can( j 
are the B candidate energy and momentum, calculated 
in the T(4S') cm. system, after performing a mass- and 
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FIG. 1: (color online). Illustration of the (pilbj cos 6^ b ) bin- 
ning. pj^ b is divided into ten bins each with a width of 0.15 
GeV/c. The cos bin width is 0.7, and fully contains the 
B + — > J/ipK + signal. The bins are centered at the kinemati- 
cally determined values for the two-body decay (blue circles). 
The small black dots in the figure are B + — > J/ipK + events 
generated by a Monte Carlo simulation. 



vertex-constrained fit to the leptons from the J / ip decay] . 

In order to determine the signal yield and charge asym- 
metry, we fit Mbc distributions after requiring \AE\ < 40 
McV. Each M bc distribution for 5.2 GeV/c 2 < M hc < 5.3 
GeV/c 2 is fitted with the sum of a Gaussian for the signal 
and an ARGUS function 14| for background. We simul- 
taneously fit the B + and B~ distributions with common 
shape parameters for the signal and background, except 
for the mean of the signal Gaussians. In the fit we assume 
that there is no asymmetry in the background; the effect 
of a possible asymmetry is included in the systematic 
error evaluation. We neglect a small contribution from 
correlated B background, which peaks at the signal posi- 
tion in the Mb c distribution. The effect of this "peaking 
background" is included in the systematic error as de- 
scribed below. Shape parameters and normalizations are 
allowed to vary in the fit. To correct for the kaon detec- 
tion asymmetry, which depends on the momentum p^ h 
and polar angle cos 8^ h of kaons in the laboratory system, 



we perform fits in ten bins of (p{^ h , cosOQ^)- The binning 
is shown in Fig. [TJ The bins are labeled 1 to 10 with in- 
creasing p^ h . We observe a total signal yield of 41,315 ± 
205 events. The M bc distribution for all bins combined 
is shown in Fig. [2] The measured raw asymmetry (A^p) 
in each bin is given in Table HJ 

We measure the kaon detection asymmetry, Af , in 
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data using Df - 
decay modes 15i i. 

u-^oi ! nx+ + n* 

The measured asymmetries of the above modes can be 
written as 



i7r+ {(j) -> K+K-) and D° - 
Here, we denote A x+ 



A Dt - A Dt 

^Ircc — ^FB 
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assuming the asymmetries are small. Here ^4fb denotes 
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FIG. 2: (color online). M bc distribution of B + J/tpK + 
candidates summed over all bins and for both B charge states 
(inset plot is on a semilog scale). The blue solid, blue dot- 
dashed, and red dashed curves are the total fit, the back- 
ground and the signal components, respectively. 



the forward-backward asymmetry due to the 7* — Z a 
interference in e + e~ — > cc and A* is the pion detec- 
tion asymmetry. We can extract — Aj? by subtract- 
ing the A® a c value from the A®° c with the assumption 

= ^4p B . This assumption is reasonable because 
the effect of the fragmentation on the forward-backward 
asymmetry is expected to be small. Possible deviations 
from this assumption are checked in data and included 
as a small contribution to the systematic error. The sub- 
traction is performed in bins of pion momentum, Pi ah , 
and polar angle in the laboratory system, cos^f!^, and 
the charmed meson's polar angle in the cm. system, 
cos#q MS (since cos6>q MS is correlated with cos0f ah and 

D + 

A F g depends on it). 

We reconstruct —> <j>ir + and D° K~ir + decays 
with charged tracks that originate from the vicinity of the 
interaction point. We require TZk > 0.6 for kaons and 
TZk < 0.4 for pions. The <j> meson candidates are selected 
with the requirement 1.00 GeV/c 2 < M K+K - < 1.04 
GeV/c 2 . To remove Df and D° mesons produced in B 
meson decays, we require the charmed meson momentum 
in the cm. system be greater than 2.5 GeV/c. The invari- 
ant mass distributions of £)+ — > <fiir + and D° — > K~tt + 
candidates after these requirements are shown in Figs. [3] 

We first obtain the A rc B c map in three-dimensional (3D) 

bins of [Pi ab , cosflj^, cos^ms)- ^ n eacn bin, ^ rc c i s ob- 
tained by fitting the reconstructed Df candidate mass 
distributions (1.89 GeV/c 2 < M 07r+ < 2.09 GeV/c 2 ). 
The fits are performed in a manner similar to those per- 
formed for B + — > J/'if)K + . D+ signals are parameter- 
ized as a sum of two Gaussian functions and a bifurcated 
Gaussian, which represents the tail of the distribution 
(for high-statistics bins only). The bin-dependent frac- 
tions of Gaussians are fixed to those obtained from Monte 
Carlo simulations. The background is parametrized as a 



first-order polynomial and its asymmetry is allowed to 

vary in the fits. Figure 2] shows the A^ e % map in bins of 
(pf ab , cos0f ab , cos0g MS ). The 3D binning (3 x 3 x 3) is 
selected to have sufficient granularity with large enough 
statistics in each bin. 

The Af values are extracted using D° — > K~ir + de- 
cays as follows: fits are performed to the D° candidate 
invariant mass distribution in 1.79 GeV/c 2 < Mx-,+ < 
1.99 GcV/c 2 with a parameterization similar to that 
used for the D+ . Here Ng c an d nD° 

are corrected 

D + 

according to the A Te s c values in bins of {Pi ah , cos8^ h , 
cos#cms)- The obtained values of N® c and N®° are 
already corrected for Ap^ and A™ + and their asymme- 
try gives — Af . The values of the kaon asymmetry Af 
are determined for different (pj^ b , cos 0^ h ) bins defined in 
Fig- HI These results are shown in Fig.[5]and the obtained 
kaon asymmetry values are quoted in the last column of 
Table |H 

Finally, the measured A™p values are corrected by 
Af + in each bin. The results arc shown in Fig. [5] 
and summarized in Table fl] We obtain Acp(B + —> 
J/ipK+) = (-0.76 ± 0.50)%, which is a weighted av- 
erage of corrected asymmetries, and in which the error 
includes only the statistical error of the A^p determi- 
nation. The weights in the averaging procedure are also 
based only on the statistical error of the A^p. Statistical 
and systematic uncertainties of A^ + are included in the 
systematic error as described below. 

TABLE I: Summary of asymmetries (in units of %). Acp 
is the corrected charge asymmetry. The final Acp values are 
weighted averages with the statistical errors from A r ^p. The 
first error in is the statistical error of the D° — s- K~tt + 

signal yield. The second error comes from the statistical errors 

D + 

in Aiei and has correlation among bins. 



Bin 


Acp 


Acp 


Af + 


1 


+0.75 


+2.30 ± 2.47 


-1.55 + 0.35 + 0.26 


2 


-1.91 


-1.04 ± 1.49 


-0.86 + 0.24 + 0.22 


3 


-2.23 


-1.65 ± 1.47 


-0.58 + 0.23 + 0.21 


4 


-0.36 


+0.16 ± 1.46 


-0.52 + 0.22 + 0.19 


5 


-0.41 


+0.01 ± 1.46 


-0.42 + 0.21 ±0.18 


6 


-2.52 


-2.19 ± 1.42 


-0.33 + 0.20 + 0.18 


7 


+1.05 


+ 1.04 ± 1.41 


+0.01 + 0.20 + 0.19 


8 


-0.14 


+0.20 ± 1.41 


-0.35 + 0.22 + 0.23 


9 


-0.23 


-0.01 + 1.57 


-0.22 + 0.24 + 0.24 


10 


-0.63 


-0.68 + 2.61 


+0.05 + 0.27 + 0.24 


Total 


-0.76 


-0.33 + 0.50 


-0.43 + 0.07 + 0.17 



Systematic errors arise from three sources: the system- 
atic uncertainty in A^p measurement, the uncertainty of 

the Af + due to Af c s c and due to A^ c . The systematic 
errors are summarized in Table [TTJ The systematic uncer- 
tainties due to the choice of binning, fit range, and mass 
windows are estimated from variations of the fit results, 
obtained by refitting the data using different choices. 
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FIG. 3: (color online). Df -> c/>7r + (a) and D° -> ^"Tr 4 " (b) 
invariant mass distribution summed over all bins including 
the charge-conjugate final states. The blue curve shows the 
results of the fit described in text, and the red dashed curve 
shows the signal. 
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FIG. 4: A^i map in bins of (p[ ab , cost9£ b , cos#cms)- The 
errors shown here are the statistical errors in the Df -> 
signal yield. The empty bins do not contain enough Df can- 

didates to obtain A re l and hence we assign the value of 0.0%. 




FIG. 5: (color online). Af cc (blue squares) and corrected A® ec 
(red circles), which corresponds to —A^ + in (p^ b , cos 0^ h ) 
bins. The errors in the corrected values of A® cc include the 
statistical errors of the A TC i . 



< 



123456789 10 
Bin# 

FIG. 6: (color online). Acp distribution in different bins. Er- 
rors are obtained from A T Qp and A^ + uncertainties summed 
in quadrature. The red dashed line shows the weighted av- 
erage of Acp and its uncertainty is shown as the light-blue 
band. The values obtained in the ten (p^ b , cos$i ab ) regions 
are consistent with each other. 



The dominant systematic error comes from the uncer- 
tainty in Arec , in which the statistical error in Df — > 
<f)ir + signal yields contributes 0.17%. The choice of 3D 
binning in (pj^ b , cos 6^, cos0q MS ) contributes 0.08%. 
The kaon detection asymmetry in <fi — > K + K~ cancels 
if the momentum distributions of K + and K~ are iden- 
tical in Df decay. We find a small difference between 
them that arises from asymmetry in the helicity angle 
distribution due to the interference between the <j> and 
S'-wave component in Df — > K + K~ir + decay. We es- 
timate the effect on our measurement to be 0.05% from 
the difference of momentum distributions in data. The 
effect of empty bins in the A le % map is estimated by set- 

D + 

ting A r ec values of empty bins to ±2% and results in 
a negligibly small contribution of 0.001%. The uncer- 
tainty in A® cc mainly comes from the statistical errors 
in the D a -> K~n+ (0.07%) and (p£ b , cos0£ b ) binning 
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(0.04%). A possible CP asymmetry in the D° -> K~ir + 
final state arises from the interference between decays 
with and without D° — D a mixing. The uncertainty is es- 
timated from the 95% confidence level upper limit on the 
CP- violating asymmetry, A® p = — y sin 8 sin 
using the world average of D° — D° mixing and CP vio- 



lation parameters [17[ and is found to be 0.01%. A^p is 



much smaller than A® p because there is no mixing be- 
tween and D~. We estimate the effect (0.01%) due to 

the possible difference between A^ and A^; we com- 
pare Afb in — > <fiir + and D + — > 4>tt + decays. In addi- 
tion to binning, the following sources are considered for 
the uncertainty in A^p . Based on the Monte Carlo sim- 
ulation, we estimate that the peaking background con- 
tributes 0.12% of the signal. The dominant contribu- 
tion comes from B + — > D°ir + with D° — > K + it~ where 
the 7r + 7r _ tracks are misidentified as £ + £~, and from 
B + -> J/ipK*(892) + . We estimate the systematic error 
(0.01%) using A CP (B+ -> D n 7r+) = 0.8 ± 0.8% and 
A CP {B+ -> J/V^*(892)+) = -4.8 ±3.3% @. The sys- 
tematic error due to possible asymmetry in non-peaking 
background (0.022%) is estimated by repeating fits al- 
lowing the background asymmetry to vary. 

The measurement and correction for kaon detection 
asymmetry is verified by repeating the whole procedure 
with different requirements on kaon identification. The 
corrected Aqp values are stable within the statistical un- 
certainty and estimated systematic error. 

Adding all systematic errors above in quadrature, the 
total systematic error is estimated to be 0.22%. 



In conclusion, using 772 x 10 6 BB meson pairs, 
we have measured the CP-violating charge asymmetry 
A C p{B+ -> J/tpK+) to be (-0.76±0.50±0.22)%, where 
the first uncertainty is statistical and second is system- 
atic. No significant evidence of CP violation is observed. 
Our measurement is consistent with the world average 



(0.9 ± 0.8)% |7[. This result significantly improves the 
precision from previous measurements at e + e~ collider 
experiments 0, [l9[ and supersedes our earlier result [2(| ■ 
The Belle result is slightly more precise than the DO mea- 
surement 0] . The two most precise measurements to date 
differ by less than two standard deviations. 
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TABLE II: 


Summary of systematic uncertainties. 


[7] 




Source 


% 


[8] 




Peaking background 


0.01 






ARGUS background 


0.02 


[9] 




Afbc bin width 


< 0.01 






(p£b,cos0ub) binning 


0.02 


[10] 


A D ' 

^rec 


Df — > d>7T Jr statistics 


0.17 




bin width 


< 0.01 


[11] 




M ~ i mass window 


0.02 






{Pil, cos0cms= cos ^iab) binning 


0.08 


[12] 




Empty bins 


< 0.01 






4> — > K K~ asymmetry 


0.05 


[13] 


A D " 


D° — > K~7v + statistics 


0.07 






M K 7r + bin width 


< 0.01 


[14] 




A / K 7 . ■ mass window 


< 0.01 






(p£b)COS0ub) binning 


0.04 


[15] 




possible Aqp 


0.01 






j4pg = Ap B assumption 


0.01 


[16] 




Total 


0.22 


[17] 
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